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Abstract.  Here  we  consider  the  results  of  an  experimental  investigation  of  the  temporal  evolution  of  plasmas  produced  by 
high  power  laser  irradiation  of  various  types  of  target  materials.  The  experiment  was  performed  at  the  LULI  Laboratory 
(Ecole  Polytechnique,  Paris).  We  have  developed  a  method  to  analyze  time-resolved  streak-camera  images  and  analyzed 
a  number  of  results  obtained  with  various  materials. 


INTRODUCTION 

The  aim  of  this  work  is  to  analyze  the  process  of  plasma  generation  and  expansion  by  high  power  laser 
irradiation  of  different  materials,  in  order  to  find  the  temporal-spatial  distribution  of  electron  plasma  density,  to  find 
the  temporal  and  spatial  profiles  of  characteristic  velocities  of  plasma  expansion  and  to  analyze  these  behaviors  vs. 
the  atomic  number  Z  and/or  the  laser  energy. 

The  investigation  of  hydrodynamics  expansion  of  laser -produced  plasmas  (in  the  intensity  range  I  ~  1014  W/cm3) 
is  of  fundamental  importance  for  several  physics  areas  such  as  inertial  confinement  fusion.  Although  several 
theoretical  models  of  plasma  expansion  were  developed  already  in  the  70’s  and  in  the  80's  [1-3]  and  many 
experiments  have  studied  this  aspect,  still  there  are  not  many  “clean”  experimental  results.  Indeed  most  experiments 
are  influenced  by  2D  effects  in  the  hydro  expansion  of  the  plasma,  which  arise  as  a  consequence  of  the  small  focal 
spots,  which  were  needed  to  produce  the  relatively  high  laser  intensity  of  interest  in  these  experiments. 

In  recent  years,  several  techniques  have  been  introduced  which  allow  the  production  of  “flat-top”  intensity 
profiles.  Hence  one  fundamental  experimental  parameter,  the  laser  intensity  on  target,  is  clearly  defined,  which  is  not 
the  case  with  the  usual  Gaussian-like  intensity  distribution  or,  even  worse,  with  focal  spots  affected  by  laser  hot 
spots.  This  allows  a  well-characterized  parametric  study  of  plasma  expansion.  Also,  2D  effects  in  plasma  expansion 
are  strongly  reduced,  producing  a  situation  much  closer  to  that  described  by  ID  theoretical  models.  One  of  such 
technique  is  that  of  Phase  Zone  Plates  (PZP  [4])  which  appeared  in  recent  years,  as  an  useful  investigation  tool,  and 
which  was  also  used  in  our  experiment. 


EXPERIMENTAL  SET  UP 

The  experimental  set  up  includes  a  Nd:glass  high  power  laser  system  with  typical  intensity  of  1014  W/cm3  and 
duration  600  ps,  a  probe  beam  (Nd:YAG  converted  to  2co)  coupled  to  an  interferometer  and  to  a  streak-camera  with 
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ps  resolution.  The  diagnostic  system  allows  the  evolution  of  the  plasma  density  profile  to  be  measured  as  a  function 
of  time.  Different  targets  (CH,  Al,  Au)  have  were  in  the  experiment,  with  different  radiative  properties.  In  particular 
plasma  expansion  with  Au  targets  should  be  strongly  influenced  by  radiation  (XUV)  transport  [5],  while,  in  the  case 
of  CH  targets,  energy  transport  should  be  determined  by  electrons  only.  Such  phenomena  are  a  key  problem  in  the 
physics  of  Inertial  Confinement  Fusion  (ICF)  driven  by  lasers. 


FIGURE  1.  Typical  experimental  set-up  used  at  LULI 


Phase  Zone  Plates  (PZP)  were  used  to  smooth  the  laser  beam  (remove  large-scale  hot  spots)  and  produce  a  flat- 
top  irradiation  profile,  reducing  2D  in  plasma  expansion  [4],  The  experimental  set-up  is  shown  in  fig.  1.  Fig.  2  shows 
typical  experimental  streak-camera  images:  a)  CH  and  b)  Al  targets  accordantly.  The  total  image  dimensions  are  375 
pm  (horizontal  scale)  and  (1.065  ns  vertical  scale). 


taiget  vacuum.  /  plasma  taiget  vacuum  /  plasma 


FIGURE  2.  Typical  interferograms  recorded  with  the  streak  camera:  a)  CH  and  b)  Al  targets.  The  total  image  dimensions  are  375 
pm  (horizontal  scale)  and  (1.065  ns  vertical  scale). 


DETERMINATION  OF  PLASMA  DENSITY 

Here  we  introduce  the  method  used  for  the  analysis  of  interferometric  images,  like  those  presented  in  fig.  2.  The 
phase  shift  in  a  point  x  of  the  interferometric  picture  from  classical  optics  is  AOa(x)  =  2nccx  /  A,  where  a  is  the  angle 
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between  the  probe  and  reference  beams  the  in  interferometer;  and  A  the  wavelength  of  the  laser  beam.  The  plasma- 


induced  phase  shift  is  A</> piasma  (x,  t)  : 
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where  Ne  is  the  density  of  free  electrons. 
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the  critical  density,  A  the  laser  wavelength.  Here  d  is  the  transversal  dimension  of  the  plasma 


(along  probe  beam,  i.e.  perpendicularly  to  x).  In  this  preliminary  analysis,  we  have  taken  d  =  400  pm,  corresponding 
to  the  FWHM  of  the  laser  focal  spot  created  by  the  PZPs.  Therefore  the  density  can  be  obtained  from  the  phase  shift 
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.In  the  limiting  case  Ne«nc,  we  get  Ne  ~  AN,.  i,  where  ANe  =  nc  A  /d,  and  where 
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i  is  the  number  of  71-phase  shifts.  Hence,  every  phase  shift  of  n  along  the  temporal  axis  corresponds  to  a  density 
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in  the  general  case,  or  ANe  =  n,  A  / d,  in  the  case 


Ne«nc. 

In  the  following  step  of  the  analysis,  we  consider  the  points  on  the  x-axis  (spatial  dimension),  which  initially,  i.e. 
before  plasma  production,  correspond  to  maxima  and  minima  of  the  interferograms,  and  we  follow  here  the  time 
evolution  of  the  density  Ne  in  steps  of  ANe. 

For  the  analysis  of  the  interferograms  we  have  applied  a  technique  of  averaging  with  time  shifts,  which  is  similar 
to  the  application  of  the  direct  Radon  transformation,  which  was  discussed  in  [6]. 


FIGURE  3  (a),  time-dependence  of  intensity  in  the  streak-camera  image  in  fig.  2a,  for  Xo= 48 1  pixels;  (b).  time-dependence 
of  intensity  in  the  streak-camera  image  in  fig.  2a,  for  x0= 48 1  after  time-shift  averaging  in  the  spatial  interval  +20  pixel.  The 
dotted  line  corresponds  to  characteristic  velocity  1.6,  the  solid  line  to  1.8.  Finally  a  Gauss-filter  by  t  as  been  applied 


ANALYSIS  OF  INTERFEROMETRICAL  IMAGES 

For  the  analysis  of  experimental  interferograms,  we  have  attempted  to  consider  a  precise  numerical  method  based 
on  averaging  with  shift  on  the  time  axis. 

Indeed,  if  consider  ione  interferogram,  as  in  fig.  2,  and  we  plot  the  intensity  corresponding  to  one  given  position 
Xo  (vertical  densitometry),  the  presence  of  a  very  strong  noise  (fig.  3  (a))  can  result  either  in  neglecting  or  shifting  of 
real  extremes  of  phase  shift,  or  in  the  appearance  of  false  extremes.  But  on  the  other  hand,  averaging  of  several  of 
such  densitometries  on  a  small  interval  Axo,  centred  around  Xo,  does  not  give  good  results,  because  it  smoothes  the 
picture. 


421 


In  alternative,  the  averaging  with  time  shifts,  with  a  characteristic  velocity  more  or  less  correspondent  to  real 
velocity,  gives  a  good  resolution  and  a  good  contrast  in  the  obtained  densitometry. 

Fig.  3  (a)  and  (b)  give  the  illustration  of  the  method  for  the  precise  determination  of  extremes  in  the 
interferometric  streak-camera  image,  shown  in  the  fig.  2  (a),  where  x0=481  pixels,  corresponding  to  the  first 
minimum  from  the  right.  The  horizontal  time  scale  is  also  in  pixels.  The  real  scale  is  512  /  512  pixels  for  375  pm  / 
1065  ps.  The  vertical  scale,  giving  the  intensity  on  the  streak-camera  image,  is  in  arbitrary  units.  How  we  can  see 
from  fig.  3  (b),  the  differences  between  locations  of  extremes  for  a  quite  wide  range  of  characteristic  time-shift 
velocities  are  not  more  than  1  pixel  (the  principle  limit),  showing  that  our  results  are  quite  reliable. 


RESULTS  AND  DISCUSSION 

Let  us  turn  back  to  fig.  2.  For  the  each  fixed  x,  we  find  3  zones:  (i)  the  first,  where  the  signal  is  only  modulated  in 
x  and  independent  from  t  (vertical  interference  fringes),  (ii)  the  second  one,  where  the  signal  significantly  depends 
on  t  (inclined  interference  fringes),  and  (iii)  third  one,  practically  uniform  where  the  interference  fringes  disappear. 
The  first,  up  to  a  time  t0(x),  correspond  to  the  situation  without  any  plasma.  At  the  time  t0,  the  first  plasma  reaches 
the  point  x  (or  at  least  a  plasma  with  a  density  sufficient  to  produce  a  measurable  phase  shift).  The  second  zone  ends 
at  the  time  tcr,  at  which  the  plasma  density  overcomes  some  critical  value,  so  that  the  probe  beam  is  refracted  away 
from  the  collecting  optics.  After  tcr,  of  course  the  fringes  disappear  and  it  becomes  impossible  to  determine  the 
plasma  density.  So,  the  zone  which  we  can  investigate  is  limited  by  the  curves  t0(x)  and  tcr(x),  which  separate  the 
region  of  non-zero  and  non-overcome  electron  density. 

Within  this  zone,  for  each  fixed  x,  we  can  further  distinguish  3  regions:  (i)  first  growth,  (ii)  constant  behaviour, 
and  (iii)  final  fast  growth.  For  a  better  analysis  of  such  behaviour,  we  consider  the  dependence  of  the  growth 
velocity  dN/dt  from  Ne  (fig.  4  shows  such  dependence  for  x=303  /mi). 
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FIGURE  4.  Velocity  growth  of  electron  density,  vs.  electron  density,  forx=303  pm,  as  deduced  from  fig.  1  a. 


Starting  from  the  typical  dependence  shown  in  fig.  4,  we  can  write: 

dNJdt  =  (ANe  /At\ ,)  f(Ne,x )  (1 ) 

where  ANe  and  Zb  corresponds  to  the  previously  introduced  density  jump  inducing  a  71-phase  shift  and  to  the 
time  difference  between  value  of  Ne  with  phase  shift  n  in  the  interferometric  picture  in  the  second  (constant)  region. 
Solving  this  differential  equation  we  can  find  the  formula  for  the  time  t(Ne,x)  ,  at  which  an  electron  density  Ne  is 
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ATm  r 

realized  in  the  point  X.  This  is  t(N  ,x)  =  tfx)  H - — 

AN.  J 


At,,  '  r  dN„ 


A Ne  J  /(Ae,x) 


where  f0  (x)  is  the  time,  where  the  first  plasma 


reaches  the  point  x,  as  before.  For  the  function  / ( Ne,x )  we  could  find  an  empirical  equation: 

A(x)  +  (1  -  A(x))  ■  (Ne  / Ncrl )  ,  Ne<  Ncrt 

f(Ne,x)  =  \  1  ,  N^KN'KN^ix)  (2) 

exp(5(x)  •  ((Ne  /  Ncr2  (x))  - 1))  ,  Ne>  Ncr2  (x) 

where  A(x)  is  an  empirical  linear  function  characterizing  of  the  growth  of  Ne  in  the  first  region;  B(x)  =  ln(l/H(x)), 
and  Ncrl  and  Ncr2  are  the  electron  densities  corresponding  to  the  beginning  and  the  end  of  the  constant  phase  in  fig. 
4.  It  is  necessary  to  remark  that  if  formally  Nrri">Nr:r2  the  function  fix)  has  no  constant  region.  We  can  finally  find  the 
dependence  of  Ne  from  t  and  x.  If  the  function  function  fix)  has  a  constant  region,  this  is: 


Ne{t,x)  =  • 


A^l  •  A(x)  (  (t  -  t0(x))  ■  (1  -  A(x))  -ANe  j  { 
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where  tcr](x)  =  t0(x)  + 


At«NcrlB(x) 


NNe  ■  (1  -  A(x)> 

while,  if  the  function  f(x)  has  no  constant  region 


tcrlix)  =  tcri(x)  +  Afii 
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where  tcr 2(x)  : 


:  t0(x)  +  At,, 


Ncrl(x) 
1  —  A(x) 


•  In  1  + 


A^cr2(x)  •  (l  ~  A(x)) 
A(x)  •  Ncrl(x) 


.  From  our  preliminary  analysis,  we  can  conclude 


that  Ncrl  is  more  or  less  independent  on  x,  Ncr2  strongly  depends  on  the  point  x.  At  is  practically  independent 
from  x  for  the  CH  target. 

From  the  calculated  values  of  Ne(t,x)  we  can  reconstruct  the  interferograms.  We  can  observe  (see  fig.  5,  which 


presents  the  comparison  between  reconstructed  interferograms  (stripes  without  noise)  and  experimental  ones 
(stripes  with  noise)  for  the  cases  presented  in  fig.  1,  CH  (left)  and  A1  (right)  targets.)  how  these  correspond  very 
well  to  the  experimental  interferograms  for  all  values  of  x  and  t. 
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FIGURE  5.  Comparison  between  reconstructed  interferograms  (stripes  without  noise)  and  experimental  ones  for  the  cases 
presented  in  fig.  1,  CH  (left)  and  A1  (right)  targets. 


CONCLUSION 

In  conclusion,  we  have  found  a  reasonable  way  to  analyze  the  interferograms  and  reconstruct  them,  and  we  have 
found  some  phenomenological  behavior.  It  should  be  stressed  that  the  values  of  Ne,  which  we  calculated,  are  NOT 
the  real  values  of  electron  density.  Indeed  they  have  been  calculated  by  assuming  a  plasma  transversal  size  d  =  400 
pm  (equal  to  the  FWHM  size  of  the  focal  spot)  and  keeping  it  constant  in  time.  Hence  in  reality  the  variation 

evidenced  in  fig.  4  rather  corresponds  to  the  evolution  of  the  line  integral  f  N e(x,y,t)dy  (y  being  the  coordinate 

along  the  probe  beam,  perpendicular  to  x)  which  we  assumed  to  be  ~(Ne  d).  A  more  proper  analysis  will  need  to 
consider  the  time  evolution  of  the  plasma  transversal  size,  through  a  2D  hydrodynamics  modeling. 
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